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Abstract 

This manual describes how to use the MatCHig code for matching the charged Higgs 
boson production processes gg —> tbH^ and gb —> A negative term, correct¬ 

ing for the double-counting between these processes, is implemented as an external 
process to Pythia, allowing the Monte Carlo generation of matched events. Re¬ 
sults from the matching were published in pQ. The code can be downloaded from 
http://www.isv.uu.se/thep/MC/matchig/. 


1 Physics motivation 

Many extensions of the Standard Model, most notably supersymmetric theories like the 
MSSM, predict the existence of a charged Higgs boson. Since there is no charged scalar 
particle in the Standard Model, the discovery of such a particle would be a clear signal 
for physics beyond the Standard Model, and would provide insight into the nature of 
the extension needed. However, in order to hnd the charged Higgs we need an accurate 
description of its production and what experimental signatures to look for. 

Charged Higgs production at hadron colliders is usually described in Monte Carlo event 
generators such as Pythia j2j and Herwig j2] using the 2 —2 process 

gb ^ tH- {gb ^ tH+) . (1) 

Here, the incoming 6-quark is described by the DGLAP evolution equations using gluon 
splitting to bb in the collinear approximation. Therefore the event also contains an outgoing 
b (6) quark. For low transverse momenta of this accompanying b quark, the 2 —2 process 
o, including initial state parton showers, describes the cross-section well, since the b quark 

density includes a resummation of potentially large logarithms of the type log , 
where /Xi? is the factorization scale. However, for large p±^b (transverse momentum of the 
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accompanying 6-quark) the collinear approximation is no longer appropriate, and one need 
instead use the exact matrix element of the 2 —> 3 process 

gg tbR- {gg tbH^) . (2) 

At low transverse momenta, the 2 —3 process can be described in terms of gluon splitting 
to bb (at order as log times the matrix element of the 2 —2 process. This means that 
we can not simply add the two processes without getting double-counting for low values 
of p±^i,. This problem is usually addressed by using only the 2 —2 process (CD when the 
accompanying b quark is not observed, and only the 2^3 process o when observing 
the accompanying b quark. However, as was shown in pp, for low transverse momenta 
{p±,b ^ 100 GeV) the latter approach underestimates the differential cross-section. There¬ 
fore, when the accompanying 6-quark is observed, it is necessary to use both the 2 —>■ 2 
and the 2 —> 3 processes together, appropriately matched to remove the double-counting. 

Detailed discussions, derivations, references and results from this matching were published 

in p. 


2 The double-counting term 


In order to remove the double-counting, we define a double-counting term ctdc; given by 
the part of the 2 —3 process (eq. (jl|)) which is already included in the 2 —>• 2 process 
(eq. (P). This term is then subtracted from the sum of the cross-sections from the two 
processes, 

O' = 02^2 + <^2^3 — O'DC • (3) 


The double-counting term is given by the leading {0{as log ^)) contribution of the 6 quark 
density to the 2 —2 process cross-section: 

da2^2 

dxidx2 


{Xi,X2) +Xi^ X2 


( 4 ) 


oac = / dxidx2 


g{Xi,flF)b'{x2,flF) 


where 6'(a:,/i|,) is the leading order 6-quark density given by 

b'{x, /i|) ~ ^ log ^ f g (5) 

with Pqg the g ^ qq splitting function, g[x,g?p) the gluon density function, pf the factor¬ 
ization scale and z the longitudinal gluon momentum fraction taken by the 6-quark. 

Including kinematic constraints due to hnite center of mass energy and hnite 6 quark mass. 
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the resulting expression for the double-counting term can be written as 


/•i dr f 

o-DC = / — / dy - —d{cose)\M2^2\ -x 

^ JUogr s ./_1 2 27r 


dzPqg (z) 


S J-i 2 

dm 

Iq2 . Q'^ + ml 

^ min ^ 


-Q: 


2 

max 


—S' (—, hi) a:25'(a^2, hi) + 

Z \ Z / 


2:2 


( 6 ) 


Here M. 2^2 is the matrix element for the 2 —>■ 2 process o, fip and fipt are the factorization 
and renormalization scales as in the 2 —> 3 process, and the kinematical variables are 
T = X 1 X 2 , Xi ^2 = s = Ts. 9 is the polar angle of the t-quark in the CM system of 

the 2 —>• 2 scattering, and Pu = s~^^y{s — mp — is the virtuality 

of the incoming 6-quark and z is identihed with the ratio of the center-of-mass energies of 
the gb (2 —2) system and the gg (2 —>• 3) system. 

The integration limits are given by 


T-min = {mt + mH±y/s 

Qlin = m - 1) - ml] - {z-^ - 1) - mlf - 4sml 

QLx = min |hl, 2 ^ 2 , 1 [s m - 1 ) - ml] + \\l[s{zm - 1) - mlf - Asml^ 

where hi 2^2 i® factorization scale of the 2 —>• 2 process, which sets the upper 
in the parton showers. 


sml, pI o^s 


(7a) 

(7b) 

(7c) 

(7d) 

limit 


3 Implementation 

The double-counting term is implemented as an external process using the Les Houches 
generic user process interface for event generators 0. This means that it can in principle 
be used together with any event generator supporting the Les Houches standard, but it 
is primarily intended for use with Pythia j2], and it uses several Pythia routines and 
common blocks. The relevant FORTRAN-hles can be downloaded from 
http://www.isv.uu.se/thep/MC/matchig/ . 

The implementation, found in the hie matchig.f, provides the subroutines UPINIT, to set 
up the external process, and UPEVNT, which supplies an event (he. a set of incoming and 
outgoing particles and momenta) and a weight for the event. Please note that since the 
double-counting contribution should be subtracted from the sum of the positive processes 
m and 0, this weight is negative for double-counting events. This means that if all three 
processes are run simultaneously in Pythia, the total cross-section will be the correctly 
matched one given by eq. ©• 
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The events from the double-counting term should likewise be subtracted in the analysis. In 
practice, this means that these events should undergo the same cuts and the same detector 
response simulation, etc., as the events from the positive contributions, and then be added 
with weight —1 to any histograms. The events from the double-counting term can be 
identihed in Pythia using the parameter construction KFPR(MSTI(1) ,2), which contains 
the external process identiher for external processes. The external process identiher for 
the double-counting process is LPRUP(l)=10000. 

Since the 2^2 process has larger cross-section than the double-counting term in all 
phase-space regions, the risk of getting a negative total cross-section in any phase-space 
point should be small given enough statistics. 

3.1 Parameters and recommended settings 

For information on how to use Pythia with external processes, please see the Pythia 
manual [2], section 9.9. An example main program, matchex.f, illustrating how to use the 
double-counting term together with Pythia, can also be found on the homepage. Please 
note that the double-counting term should be run together with the internal processes 
gb tH^ (ISUB=161) and gg —^ tbH^ (ISUB=401) to get matched events. (Note that the 
latter process (ISUB=401) was implemented in Pythia from version 6.223.) 

The following Pythia parameters are used in matchig.f; 


PMAS (in the common block pydat2) is used to get masses and widths for the 

charged Higgs boson (PMAS (37,1), PMAS(37,2)),top quark (PMAS (6,1), 
PMAS (6, 2)), bottom quark (PMAS (5,1)) and boson (PMAS(24,1)). 

PARU(141) (in PYDATl) is used to get the value of tan/3. 

PARU(102) (in PYDATl) is used to get the value of sin^ 9w- 

MSTP(39) (in PYPARs) is used to choose factorization and renormalization scales as 

for the 2 —3 process. The following choices are possible: MSTP(39) = 

3, 5, 6 or 8. The options MSTP(39) = 6 and 8 are new from Pythia 
6.226, please see [S] for details. Our recommended setting is MSTP(39) = 
8, corresponding to MSTP(32) = 12. Note that in this case PARP(193) and 
PARP(194) must also be set, see below. 

MSTP(32) (in PYPARs) is used to choose the maximum scale for parton showers, i.e. 

the maximum transverse momentum of the accompanying 6-quark as for 
the 2 —»■ 2 process (/if, 2^2 in eq. ((Zj)). The following choices are possible: 
MSTP(32) =4, 11 or 12. The options MSTP(32) =11 and 12 are new 
from Pythia 6.226, please see |3] for details. Our recommended setting is 
MSTP(32) = 12, corresponding to MSTP(39) = 8. Note that in this case 
PARP(193) and PARP(194) must also be set, see below. 
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PARP(34) (in PYPARs) is used to modify the factorization and renormalization scales 
as in Pythia. Note however that it does not affect the maximum trans¬ 
verse momentum of the accompanying 6-quark. 

PARP(67) (in PYPARs) is used to modify the maximum scale for parton showers, i.e. 
the maximum transverse momentum of the accompanying 6-quark 
in eq. dlD), as in Pythia. 

PARP (48) (in PYPARs) is used to cut the tails of the Breit-Wigner distribution for the 
charged Higgs mass at masses outside mH± ±PARP(48)P h±, as in Pythia. 

The values of Oem, the parton distributions and the running top and bottom quark 
masses are calculated using the corresponding Pythia subroutines. The top quark and 
charged Higgs masses are Breit-Wigner distributed with mass-dependent widths in the 
same way as in the internal gb —> tH^ and gg —> tbH^ processes, using the Pythia 
subroutine PYWIDT. Suppression of the cross-section due to closed decay channels of the 
charged Higgs boson and top quark is done as in the Pythia processes. 

Factorization scale: As shown in P, there are strong arguments that a non-standard 
factorization scale gip ~ should be used in charged Higgs production. The 

option to have hxed factorization and renormalization scales of choice was added to Pythia 
from version 6.226 jH]. This is done by setting MSTP(32)=12 (for 2 —2 processes) and 
MSTP(39)=8 (for 2 —3 processes). Then one needs to set PARP(193) to the square of 
the factorization scale (/i^), and PARP (194) to the square of the renormalization scale 
(usually taken to be 

On-shell parton showers: In order to be able to compare the differential cross-section for 
the 2 —2 process (HD, the 2 —>• 3 process m and the double-counting term on parton level, 
one should use on-shell initial-state parton showers, set by MSTP(63)=0. In a realistic case, 
using jet-hnding algorithms, there should however be no observable difference compared 
to having the parton showers off-shell. 


4 Final comments 

The code described in this document was written for research conducted together with 
Johan Rathsman. If you use this code, please cite PJ. 
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